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We describe the measurement of the depth of maximum, J*L max , of the longitudinal development 
of air showers induced by cosmic rays. Almost four thousand events above 10 18 eV observed by 
the fluorescence detector of the Pierre Auger Observatory in coincidence with at least one surface 
detector station are selected for the analysis. The average shower maximum was found to evolve with 
energy at a rate of (106^) g/cm 2 /decade below io 18 ' 24±0 05 e V and (24±3) g/cm 2 /decade above 
this energy. The measured shower-to-shower fluctuations decrease from about 55 to 26 g/cm 2 . The 
interpretation of these results in terms of the cosmic ray mass composition is briefly discussed. 

PACS numbers: 96.50.sd,13.85.Tp,98.70.Sa 



Introduction - The energy dependence of the mass 
composition of cosmic rays is, along with the flux and ar- 
rival direction distribution, an important parameter for 
the understanding of the sources and propagation of cos- 
mic rays at very high energy. There are several models 
that describe the observed flux of cosmic rays very well, 
but each of these models has different assumptions about 
the cosmic ray sources and correspondingly predicts a 
different mass composition at Earth. For example, the 
hardening of the cosmic ray energy spectrum at ener- 
gies between 10 18 eV and 10 19 eV, known as the 'ankle', 
is presumed to be either a signature of the transition 
from galactic to extragalactic cosmic rays or a distor- 
tion of a proton-dominated extragalactic spectrum due 
to energy losses [H. Moreover, composition information 
may eventually help to decide whether the flux suppres- 
sion observed above 4-10 19 eV Q is due mainly to the in- 
teraction of cosmic rays with the microwave background 
or a signature of the maximum injection energy of the 
sources 0. 

Due to the low flux at these energies, the composition 
of cosmic rays cannot be measured directly, but has to 
be inferred from observations of extensive air showers. 
The atmospheric depth, X max , at which the longitudinal 
development of a shower reaches its maximum in terms 
of the number of secondary particles is correlated with 
the mass of the incident cosmic ray particle. With the 
generalization of Heitler's model of electron-photon cas- 
cades to hadron-induced showers and the superposition 
assumption for nuclear primaries of mass A, the average 
depth of the shower maximum, ( X max ) , at a given energy 
E is expected to follow Q 



The change of (A" max ) per decade of energy is called 
elongation rate Q, 



(X max )=a(\nE-(\nA))+(3,) 



(1) 



where (In A) is the average of the logarithm of the pri- 
mary masses. The coefficients a and (3 depend on the 
nature of hadronic interactions, most notably on the mul- 
tiplicity, elasticity and cross-section in ultra-high energy 
collisions of hadrons with air, see e.g. [5]. Although 
Eq. ([TJ is based on a simplified description of air showers, 
it gives a good description of air shower simulations with 
energy-independent parameters a and j3 in the energy 
range considered here, see @. Only physics processes 
not accounted for in currently available interaction mod- 
els could lead to a significant energy dependence of these 
parameters. 



D 



10 



d(X n 



dlgE 



a 1 



d(\nA) 
dlnE 



ln(10), (2) 



and it is sensitive to changes in composition with en- 
ergy. A complementary composition-dependent observ- 
able is the magnitude of the shower-to-shower fluctua- 
tions of the depth of maximum, RMS(AT max ), which is 
expected to decrease with the number of primary nucle- 
ons A (though not as fast as 1/y/A [8() and to increase 
with the interaction length of the primary particle. 

At ultra high energies, the shower maximum can be 
observed directly with fluorescence detectors. Previously 
published AT max measurements 0, [To| focused mainly on 
(I mra ) as a function of energy and had only limited 
statistics above 10 19 eV. 

Here we present a measurement of both (X max ) and 
RMS(A" max ) using high quality and high statistics data 
collected with the southern site of the Pierre Auger Ob- 
servatory [lH . The Observatory is located in the province 
of Mendoza, Argentina and consists of two detectors. 
The surface detector (SD) array comprises 1600 water- 
Cherenkov detectors arranged on a triangular grid with 
1500 m spacing that cover an area of over 3000 km 2 . The 
water-Cherenkov detectors are sensitive to the air shower 
components at ground level. The fluorescence detector 
(FD) consists of 24 optical telescopes overlooking the ar- 
ray, which can observe the longitudinal shower develop- 
ment by detecting the fluorescence and Cherenkov light 
produced by charged particles along the shower trajec- 
tory in the atmosphere. 

Data Analysis. - This work is based on air shower data 
recorded between December 2004 and March 2009. Only 
events detected in hybrid mode jl2j are considered, i.e. 
the shower development must have been measured by the 
FD, and at least one coincident SD station is required to 
provide a ground-level time. Using the time constraint 
from the SD, the shower geometry can be determined 
with an angular uncertainty of 0.6° [l3j]. The longitu- 
dinal profile of the energy deposit is reconstructed [14 1 
from the light recorded by the FD using the fluorescence 
and Cherenkov yields and lateral distributions from [lj| . 
With the help of data from atmospheric monitoring de- 
vices [HI the light collected by the telescopes is corrected 
for the attenuation between the shower and the detector 
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AX max /V2 [g/cm 2 ] 

FIG. 1: Difference between X max measured in showers simul- 
taneously at two FD stations ({lg(_B/eV)> = 19.1). The X max 
resolution is displayed as a function of energy in the inset. 

and the longitudinal shower profile is reconstructed as 
a function of atmospheric depth. -V max is determined 
by fitting the reconstructed longitudinal profile with a 
Gaisser-Hillas function [17] . 

An unbiased set of high quality events is selected with 
the statistical uncertainty of the reconstructed X max be- 
ing comparable to the size of the fluctuations expected 
for nuclei as heavy as iron (as 20 g/cm 2 ) and small sys- 
tematic uncertainties as explained in the following. 

The impact of varying atmospheric conditions on the 
A^ max measurement is minimized by rejecting time peri- 
ods with cloud coverage and by requiring reliable mea- 
surements of the vertical optical depth of aerosols. Pro- 
files that are distorted by residual cloud contamination 
are rejected by a loose cut on the quality of the profile 
fit (x 2 /Ndf<2.5). We take into account events only with 
energies above 10 18 eV where the probability for at least 
one triggered SD station is 100%, irrespective of the mass 
of the primary particle [l8| ■ The geometrical reconstruc- 
tion of showers with a large apparent angular speed of the 
image in the telescope is susceptible to uncertainties in 
the time synchronization between FD and SD. Therefore, 
events with a light emission angle towards the FD that 
is smaller than 20° are rejected. This cut also removes 
events with a large fraction of Cherenkov light. The en- 
ergy and shower maximum can be reliably measured only 
if X max is in the field of view (FOV) of the telescopes 
(covering 1.5° to 30° in elevation). Events for which only 
the rising or falling edge of the profile is detected are 
not used. Moreover, we calculate the expected statisti- 
cal uncertainty of the reconstruction of X max for each 
event, based on the shower geometry and atmospheric 
conditions, and require it to be better than 40 g/cm 2 . 

The latter two selection criteria may cause a selection 
bias due to a systematic undersampling of the tails of 
the true X max distribution, since showers developing very 




E[eV] 



FIG. 2: (X max ) as a function of energy. Lines denote a fit 
with a broken line in IgE. The systematic uncertainties of 
(Xmax) are indicated by a dashed line. The number of events 
in each energy bin is displayed below the data points. HiRes 
data [Tot ] are shown for comparison. 

deep or shallow in the atmosphere might be rejected from 
the data sample. To avoid such a bias in the measured 
(X max ) and RMS(X roax ) we apply fiducial volume cuts 
based on the shower geometry that ensure that the view- 
able X max range for each shower is large enough to ac- 
commodate the full X may r distribution [19j. 

After all cuts, 3754 events are selected for the X max 
analysis. The X max resolution as a function of energy 
for these events is estimated using a detailed simulation 
of the FD and the atmosphere. As shown in the inset 
of Fig. [TJ the resolution is at the 20 g/cm 2 level above a 
few EeV. The difference between the reconstructed A" max 
values in events that had a sufficiently high energy to 
be detected independently by two or more FD stations 
is used to cross-check these findings. As can be seen in 
Fig. [TJ the simulations reproduce the data well. 

Results and Discussion. - The measured (A max ) and 
RMS(A" max ) values are shown in Figs. [2] and [3] We use 
bins of AlgE = 0.1 below 10 EeV and A lg E = 0.2 above 
that energy. The last bin starts at 10 194 eV, integrating 
up to the highest energy event (E = (59 ± 8) EeV). The 
systematic uncertainty of the FD energy scale is 22% [l8| . 
Uncertainties of the calibration, atmospheric conditions, 
reconstruction and event selection give rise to a system- 
atic uncertainty of < 13 g/cm 2 for (V max ) and <6g/cm 2 
for the RMS. The results were found to be independent 
of zenith angle, time periods and FD stations within the 
experimental uncertainties. 

A fit of the measured (A" max ) values with a con- 
stant elongation rate does not describe our data 
(x 2 /Ndf=34.9/ll), but as can be seen in Fig. [2 us- 
ing two slopes yields a satisfactory fit (x 2 /Ndf=9.7/9) 
with an elongation rate of (106^oi) g/cm 2 /decade below 
10 is.24±o.05 eV and (24±3) g/cm 2 /decade above this en- 
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FIG. 3: (X m ax) and RMS(X max ) compared with air shower simulations [20] using different hadronic interaction models [21 



ergy. If the properties of hadronic interactions do not 
change significantly over less than two orders of magni- 
tude in primary energy (< factor 10 in center of mass 
energy), this change of AD W =(82+ 21 ) g/cm 2 /decade 
would imply a change in the energy dependence of the 
composition around the ankle, supporting the hypothe- 
sis of a transition from galactic to extragalactic cosmic 
rays in this region. 

The (X max ) result of this analysis is compared to the 
HiRes data [lCf in Fig. [2] Both data-sets agree well 
within the quoted systematic uncertainties. The x 2 /Ndf 
of the HiRes data with respect to the broken-line fit de- 
scribed above is 20.5/14. This value reduces to 16.8/14 
if a relative energy shift of 15% is applied, such as sug- 
gested by a comparison of the Auger and HiRes energy 
spectra [2j. 

The shower-to-shower fluctuations, RMS(A max ), are 
obtained by subtracting the detector resolution in 
quadrature from the width of the observed -X max dis- 
tributions resulting in a correction of <6g/cm 2 . As can 
be seen in the right panel of Fig. [31 we observe a de- 
crease in the fluctuations with energy from about 55 to 
26g/cm 2 as the energy increases. Assuming again that 
the hadronic interaction properties do not change much 
within the observed energy range, these decreasing fluc- 
tuations are an independent signature of an increasing 
average mass of the primary particles. 

For the interpretation of the absolute values of (A max ) 
and RMS (.Xmax) a comparison to air shower simulations 
is needed. As can be seen in Fig. [31 there are considerable 
differences between the results of calculations using dif- 
ferent hadronic interaction models. These differences are 
not necessarily exhaustive, since the hadronic interaction 
models do not cover the full range of possible extrapola- 
tions of low energy accelerator data. If, however, these 
models provide a realistic description of hadronic inter- 
actions at ultra high energies, the comparison of the data 
and simulations leads to the same conclusions as above, 



namely a gradual increase of the average mass of cosmic 
rays with energy up to 59 EeV. 
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